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The radiative return has proven to be a competitive method for the precise measurement of the hadronic cross 
section, detailed studies of hadronic interactions, and even discoveries of new resonances. The most recent and 
future devolopments of the Monte Carlo event generator PHOKHARA are highlighted, and the impact of the 
radiative return measurements on the (g — 2) M puzzle is discussed. 



1. Introduction 

Electron-positron annihilation into hadrons is 
one of the basic reactions of particle physics, rel- 
evant for the understanding of hadronic interac- 
tions. The low energy region is crucial for pre- 
dictions of the hadronic contributions to a M = 
(g — 2)^/2, the anomalous magnetic moment of 
the muon, and to the running of the electromag- 
netic coupling from its value at low energy up to 
M z . 

The traditional way of measuring the hadronic 
cross section, the energy scan, needs dedicated ex- 
periments. An alternative and advantageous way 
is the radiative return method. This method al- 
lows for a simultaneous extraction of the hadronic 
cross section over a wide energy range in an ho- 
mogeneous data set, profiting from the data of all 
high luminosity meson factories. 

The radiative return method relies on the ob- 
servation that the cross section of the reaction 
e + e~ — ► hadrons+photons, with photons emit- 
ted from the initial lcptons (ISR), see Fig. 1(a), 
factorizes into a function H, fully calculable 
within QED, and the cross section of the re- 
action e + e~ — > hadrons. Thus from the mea- 
sured differential cross section of the reaction 
e + e~ — > hadrons + photons as a function of 
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the hadronic invariant mass one can evaluate 
cr(e + e _ — > hadrons) once the radiator function 
H is known. The radiative return method al- 
lows for the extraction of the hadronic cross sec- 
tion from the production energy threshold of a 
given hadronic channel almost to the nominal 
energy of the experiment. The smaller cross 
section of the radiative process as compared to 
the process without photon emission has to be 
compensated by higher luminosities. That re- 
quirement is met by meson factories (DAPHNE, 
CLEO, BABAR, BELLE). All of them were built 
for other purposes than the measurement of the 
hadronic cross section, but their huge luminosities 
provide with data samples large enough for very 
accurate measurements of interesting hadronic 
channels and/or give information on rare chan- 
nels, which were not accessible in scan experi- 
ments. 

Two representative examples of such measure- 
ments are the very accurate pion form factor 
extraction by the KLOE Collaboration [1], and 
the measurement of a(e + e~ — > 3 pions) by the 
BABAR Collaboration [2], where it was shown 
that the old DM2 scan data were too low at high 
values of the invariant mass of the hadronic sys- 
tem. For a review of BABAR results see Ref. [3] . 

2. The PHOKHARA event generator 

In realistic experimental situations, where 
sophisticated event selections are used, one 
needs a Monte Carlo event generator of 
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Figure 1. Leading order contributions to the reac- 
tion e + e~ — > 7r + 7r _ 7 from ISR (a) and FSR (b). 



the measured process. To meet that re- 
quirement the PHOKHARA event genera- 
tor (http://ccrn.ch/german.rodrigo/phokhara/ 2 ) 
was constructed. PHOKHARA started from the 
EVA generator [4], where the structure function 
method was used to model multi-photon emis- 
sion. The physical accuracy of EVA was how- 
ever far from the demanding experimental accu- 
racy for the measurement of the pion form fac- 
tor. PHOKHARA relies instead on exact matrix 
elements at next-to- leading order (NLO), namely 
it includes one loop radiative corrections to one- 
photon radiation and emission of two real hard 
photons. The accuracy of the simulation has been 
estimated to be of the order of 5 per mil from 
ISR [5]. 

The first version of PHOKHARA [5] was de- 
signed to run with tagged photon configurations. 
Radiative corrections necessary for photon emis- 
sion at small angles were calculated afterwords [6] 
and implemented into the event generator [7]. 
The important issue of final state emission was 
addressed subsequently [8]. In parallel the gen- 
erator was being extended to allow for the gen- 
eration of more hadronic channels. The present 
version of the program simulates the production 
of not only a pair of pions or muons, but includes 
also the simulation of K + K~ , K°K°, events with 
three [9] and four pions [7] , and nucleon pairs pp 
and nn [10] . 

All that allowed for building the state-of-thc 
art event generator. The proper implementation 

2 Mirror link at http://ific.uv.es/~rodrigo/phokhara/ 



of the radiative corrections as well as the hadronic 
currents is guarantied by extensive tests. Com- 
parisons with the KKMC [11] Monte Carlo event 
generator have been performed [12] leading to an 
excellent agreement. The comparison is however 
limited to muons in the final state. Higher order 
effects, that can be seen as a difference between 
exponentiated and non-exponentiated matrix el- 
ements reach at most 2 per mile with the excep- 
tion of the region where the hadronic system has 
an invariant mass very close to the nominal en- 
ergy of the experiment. There, soft multi-photon 
emission plays an important role and thus expo- 
nentiation is necessary. This region is however 
out of the region of interest for radiative return 
measurements [5], 

The impresive amount of new data provided in 
particular by B-factories requires further efforts 
to improuve the accuracy of the event generator, 
as well as to implement new hadronics channels. 
The latter requires a fairly good paramctrization 
of various form factors. 

3. Final state radiation and radiative <fi- 
decays 

Final state radiation (FSR), see Fig. 1(b), is 
the main background for radiative return mea- 
surements. The situation at S-factories is how- 
ever completely different from the one at the Ol- 
factory DAPHNE. In the former case the region 
of hadronic masses below 4 GeV, which is of the 
utmost physical interests, lays far from the nom- 
inal energy of the experiments. Thus an emission 
of a very hard photon is required to reach it. As 
a result the typical kinematic configuration of an 
event consists of a photon emitted back-to-back 
to the hadronic system. That provides a natu- 
ral suppression of FSR contributions, which are 
large for photons emitted parallel to the direction 
of the charged hadrons in the final state. 

At the (^-factory the physically interesting re- 
gion is not so far from the nominal energy of 
the experiment, and that natural suppression of 
FSR do not hold. Strategies should be stab- 
lished to either suppress FSR through kinemat- 
ical cuts, or to control the uncertainty due to 
the model dependence of the simulation. The 
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KLOE analysis [1,13] of the pion form factor from 
events with untagged photons emitted at small 
polar angles follows the first strategy. When 
photons in the forward-backward directions and 
pions in the central region are selected FSR is 
easily reduced to less than 1%. The price to 
pay however is that the region close to thresh- 
old, M n7V < 590 MeV, is also suppressed as pions 
are produced in this kinematical region essentially 
back-to-back to the ISR photon and therefore at 
very small polar angles that scape from detec- 
tion. For the case of untagged photons a specific 
background, e+e~ — ► 7r + 7r~e + e~ , has to be also 
taken into account [14,15] as the final leptons are 
not vetoed. 

KLOE is now perfoming a complementary 
analysis with photons tagged at large polar an- 
gles [13] that will provide accurate data in the 
threshold region. Due to the factor 1/s 2 of the 
dispersion integral for a M , that low energy region 
is highly enhanced; contributing to around 20% 
to the total integral. Therefore its relevance. 

At large photon polar angles FSR and cf) de- 
cays: <j) — > 7T + 7r~7r° and — ► /07 — ► 7r + 7r~ 7 [16]; 
become important. The background channel <f> — ► 
7T + 7r~7r can be eliminated through dedicated se- 
lection cuts, but FSR and <f> — > 7r + 7r~7 have to 
be subtractred relying on Monte Carlo. Another 
possibility to eliminate <fi decays is to run off- 
peak. DAPHNE has taken data off-resonance at a 
center-of-mass energy of = 1 GeV that will al- 
low to reduce the systematic errors in the thresh- 
old region [13]. 

The main tool to test the model dependence 
of photon emission from the final state pions and 
radiative 0-decays is the charge asymmetry. For 
events with emission of one real photon the two- 
pion state is produced with charge conjugation 
C = — 1 and odd orbital angular momentum 
when the real photon is emitted from the initial 
state, and with C = +1 and even orbital angu- 
lar momentum when the real photon is emitted 
from the final state. As a result, the ISR-FSR 
interference is odd under the exchange ir + <-> n~ 
and integrates to zero for charge blind event se- 
lections. At the same time it is the only source of 
the charge asymmetry and as such allow to test 
the FSR model. 



As shown in [17], the charge asymmetry has 
large analyzing power and can provide informa- 
tion allowing for distinguishing between differ- 
ent models of the radiative <j) ^^1 decay [18]. 
PHOKHARA have adopted two models describ- 
ing the decays cf> — > ir + ir~j and <ft ~ * 

7r 7r° 7 : 

the "no structure" model [19] and the K + K~ 
model [20]. Again by appropriate event selec- 
tions one can suppress those contributions or en- 
hance them as for other sources of FSR emission. 
Other contributions that might be important in 
the threshold region beyond the model currently 
used in PHOKHARA to describe FSR (sQED + 
vector dominance + radiative (f> decays) have been 
advocated in Ref. [21]. 

The reaction e + e~ — ► Tr + n~-f, with the photon 
emitted from the pions, does contribute also to 
dispersion integrals for the evaluation of a M and 
olqed- In the former case its theoretically esti- 
mated value [8] is of the size of the theoretical 
uncertainty and thus numerically important. As 
its theoretical estimations are not reliable it has 
to be measured. PHOKHARA also includes the 
simulation of events with simultaneous emission 
of one photon from the initial state and another 
from the final state. 

4. Three-pion channel 

The channel e + e~ — > 7r + 7r~7r° has been re- 
cently implemented in PHOKHARA [9]. The 
model for the form factor is based on generalized 
vector dominance with isospin I = 0, 1 compo- 
nents: 

jem,3« = ( 7r +( 9+ ) 7r - (g+)7r O (9o) | jem | Q) 

= W79+<£<7o Fa^il+'V-'Vo) • (1) 

A global fit has been performed with contribu- 
tions from w(782), J = w(1420), u" = w(1650), 
0(1020), p(770), p' = p(1450) and p" = p(1700) 
resonances. The model provides a very good de- 
scription of the total cross section and also pre- 
dicts in good agreement with experiment the de- 
cay width T(ir° — > 77), the slope parameter of 
7T° — > 77* , and the radiative vector meson decays 
p — > 7r°7 and cf) — ► 7r°7, but is however in con- 
flict with u) — > 7r° 7 . There is still room from im- 
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provements once information on subdistributions 
is included in the fit. 

5. Nucleon form factors 

Another example of the power of the radiative 
return method is the measurement of the mag- 
netic and the electric nucleon form factors in the 
time-like region. This measurement was first pro- 
posed in Ref. [10]. The recent measurement by 
the BABAR Collaboration [22] of the magnetic 
and electric proton form factor shows a clear evi- 
dence for a ratio \Ge/Gm\ > 1 just above thresh- 
old with respect to to previous analysis. It is also 
interesting to note that the measurement of the 
relative fase between Ge and Gm requires access 
to the nucleon spin [10] . 

6. Production and decay of J/tp and other 
narrow resonances 

J/ip resonances are copiously produced at B- 
factories. There is a strong demand from experi- 
mental groups to include this kind a new channels 
in our program, to analyse its reach phenomenol- 
ogy with high precision. Other narrow resonances 
are of interest as well. 

The BABAR Collaboration [23] has reported 
the discovery of a new state, the Y(4260) res- 
onance, by using radiative return events in the 
Tr + TT~J/ip channel. These results have been con- 
firmed by CLEO using energy-scan [24], and 
ISR [25] data, as well as by the BELLE Collabo- 
ration [26]. 

7. One-loop corrections to muon produc- 
tion and higher order radiative correc- 
tions 

When the first version of PHOKHARA was 
constructed the two-pion channel was the most 
interesting as it enters the prediction for the 
anomalous magnetic moment of the muon. By 
that time ISR radiative corrections where intro- 
duced in PHOKHARA through a leptonic tensor: 

^ISR = a2 [ fl 00 9^ + a ll P\P\ + a 22 P2P2 

(2) 



where pi (p 2 ) are the four-momentum of the in- 
coming positron (electron), and scalar 
coefficients where the antisymmetric imaginary 
contribution proportional to a_i appears first at 
NLO. While ISR corrections are independent of 
the hadronic channel, FSR and the ISR-FSR in- 
terference do depend and have to be calculated 
for each channel independently. 

When more and more hadronic channels are of 
interest it seems more convenient to introduce ra- 
diative corrections at the amplitude level by using 
the helicity amplitude formalism. Interferences 
between different amplitudes are then obtained 
automatically without further analytical compu- 
tations. The radiative one-loop amplitude can be 
factorized into three components: 

\A) = \A) 1SR + \A) FSR + |A) 2r , (3) 

where the last one steams from the exchange of 
two virtual photons between the initial and the 
final state. 

Muon pair production is not only important 
for the normalization of the i?-ratio, but being a 
very clean process that can be calculated in QED 
can be used for luminosity monitoring at e + e~ 
machines. Therefore the importance of having 
very accurate predictions. We are now complet- 
ing the calculation of ISR radiative corrections at 
the amplitude level and FSR for the muon chan- 
nel. The two-photon exchange amplitude will be 
calculated subsequently. Note that the latter is 
not even known for Bhabha scattering. 

The estimated accuracy of PHOKHARA from 
ISR is of the order of 5 per mil. Although this 
is a very conservative estimate a better accuracy 
might be needed for future experiments; higher 
luminosity B-factories or even the International 
Linear Collider (ILC). Two loop corrections, at 
least in the leading log approximation, will reduce 
this uncertainty to at least 1 to 2 per mil. 

8. Interplay between e+e and tau data 

The hadronic vacuum polarization contribution 
to the SM prediction of the anomalous magnetic 
moment of the muon is obtained though a disper- 
sion integral over the e+e~ hadronic cross sec- 
tion (once ISR and vacuum polarization correc- 
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Table 1 

Contribution of the pion form factor to aJj ad ' LO 
(in units of 1CP 11 ) from the different energy re- 
gions. Numbers from Ref. [33] (* my own esti- 
mate). 



390-520 MeV 




SND 


478.0 ±17.3 ±6.9(18.6) 


CMD-2 


461.7 ±9.8 ±3.2(10.3) 


520-600 MeV 




SND 


425.0(27)* 


600-960 MeV 




SND 


3768 ± 13 ± 47(48) 


CMD-2 


3771 ± 19 ± 27(33) 


KLOE 


3756±8±49(50) 



tions are subtracted) normalized to the point-like 
muonic cross section a = 47ra 2 /3s, where the 
low energy data has the largest weight. This re- 
quires a complicated unfolding of radiative correc- 
tions. Alternatively, conservation of vector cur- 
rents (CVC) and isospin symmetry allows to use 
the hadronic vector spectral function of tau lep- 
ton decays as input to this dispersion integral. At 
present, the SM prediction obtained from tau and 
e + e~ data do not agree to each other. This might 
be an indication of unaccounted isospin breaking 
contributions [27] but it is at the same time rather 
puzzling because precisely the tau based analysis 
agrees with the measurement of the anomalous 
magnetic moment from the BNL experiment [28], 
while the e+e~ based analysis doesn't. This dis- 
crepancy is of the order of three standard de- 
viations [29]. Recent preliminary BELLE [30] 
tau data seem to point however towards a bet- 
ter agreement between e+e~ and tau data. The 
branching ratio of r — > v T ir + ir° calculated from 
e+e~ data neither agrees [29]. 

The situation among e+e~ data is also puz- 
zling. The integral over the energy scan 
(SND [31] and CMD-2 [32]) and the radiative re- 
turn (KLOE [1]) data, and therefore the predic- 
tion for a^, agrees within 0.5er. This agreement 
is due to a tricky compensation of discrepancies 
in the shape of the energy distribution above and 
below the p-peak. 



A careful regard to these data tell us that the 
three experiments do overlap indeed only in the 
energy region above 600 MeV, while in the thresh- 
old region only data from the Novosibirsk exper- 
iments are available. As expected the statistical 
error of radiative return data is much better than 
of energy scan. CMD-2 has the better systematic 
error (0.6%), while the total systematic error of 
SND and KLOE are comparable (1.3%). The new 
KLOE analysis [13] will however reduce it to less 
than 1%. While SND and CMD-2 experiments 
agree to each other in the full energy range, they 
can not be considered completely independent as 
both analysis use the same radiative correction 
package. Furthermore, the new CMD-2 data do 
not cover the region between 520 and 600 MeV. 
This makes the ongoing KLOE large photon po- 
lar angle analysis particularly interesting because 
the agreement of the integral is very unlikely to 
happen in the threshold region. Excluding a set 
of data from the prediction of a M , until these dis- 
crepancies are solved, might lead to a biased re- 
sult. 

9. Summary 

The radiative return has proven to be a com- 
petitive method for the precise measurement of 
the hadronic cross section, detailed studies of 
hadronic interactions, and even discoveries of new 
resonances. Many interesting problems, for ex- 
ample a proper modeling of the hadronic current 
of multi-meson final states, FSR simulation for 
more than two-pions, modeling of narrow reso- 
nances and many others not mentioned await still 
for detailed theoretical investigations. 

New data from the KLOE experiment at small 
polar photon angles with a total systematic un- 
certainty below 1% and in particular in the 
threshold region from the large photon polar an- 
gle analysis [13], the long awaited pion form factor 
measurement at _B-factories, but also future new 
data from the energy scan [33] should help to clar- 
ify the discrepancies between different e+e - data 
sets. To a great extent independent test of ra- 
diative corrections are also needed. One should 
remember that even the direct measurement of 
the ratio oVtt/oVv, where most systematics are 
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expected to cancel, requires a careful treatment 
of radiative corrections and unfolding, due to the 
fact that radiative corrections affect differently 
the angular distributions of pions and muons [5] , 
and because the i?-ratio entering the dispersion 
integral for a M is normalized to the point-like cross 
section of the muon but not to the physical cross 
section. All together is necessary to confirm the 
discrepancy of the (g — 2) M measurement with re- 
spect to the Standard Model prediction. Under- 
standing tau data, whether the agreement with 
(g — 2)^ is by accident or not, is also crucial for 
this purpose. 
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